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ABSTRACT 

Proton-rich material in a state of nuclear statistical equilibrium (NSE) is one of the least studied 
regimes of nucleosynthesis. One reason for this is that after hydrogen burning, stellar evolution 
proceeds at conditions of equal number of neutrons and protons or at a slight degree of neutron- 
richness. Proton-rich nucleosynthesis in stars tends to occur only when hydrogen-rich material that 
accretes onto a white dwarf of neutron star explodes, or when neutrino interactions in the winds from 
a nascent proto-neutron star or collapsar-disk drive the matter proton-rich prior to or during the 
nucleosynthesis. In this paper we solve the NSE equations for a range of proton-rich thermodynamic 
conditions. We show that cold proton-rich NSE is qualitatively different from neutron-rich NSE. 
Instead of being dominated by the Fe-peak nuclei with the largest binding energy per nucleon that 
have a proton to nucleon ratio close to the prescribed electron fraction, NSE for proton-rich material 
near freeze-out temperature is mainly composed of 56 Ni and free protons. Previous results of nuclear 
reaction network calculations rely on this non-intuitive high proton abundance, which this paper will 
explain. We show how the differences and especially the large fraction of free protons arises from the 
minimization of the free energy as a result of a delicate competition between the entropy and nuclear 
binding energy. 

Subject headings: nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

Recently, a new nucleosynthesis process was invented 
to explain the produ ction of the proton- rich isotopes 
92 ' 94 Mo and 96 ' 98 Ru (jFrohlich et al.ll200l . In this so 
called i/p-process, matter ejected from the surface layers 
of a proto-neutron star is exposed to a strong neutrino 
flux which results in the following weak interactions: 
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The mass difference between the neutron and the pro- 
ton causes the neutrino interactions to be dominated by 
the forward reaction of eq. ([1]) , which drives the nuclear 
matter proton-rich. The matter then expands and as- 
sembles into NSE, which is not only proton-rich (i.e. 
Y e > 0.5) but actually contains a large mass fraction 
of free protons. The forward reaction of eq. @ con- 
verts some of the free protons to neutrons, which al- 
lows nuclear matter to move past the "bottle neck" nu- 
cleus 64 Ge via the fast 64 Ge(n,p) 64 Ga reaction and the 
aforementioned proton-rich isotopes are synthesized dur- 
ing the freeze out phase from the proton-rich NSE state 
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culate nucleosynthesis in similar environments. 

We show that the NSE mass fractions exhibit a great 
degree of symmetry across the line Y e = 0.5 for high tem- 
peratures where the composition is dominated by free 
nucleons and 4 He. For colder temperatures, (Tg < 6.0 
and p ~ 10 7 g cm -3 ), there are hardly any free neu- 
trons present in the NSE state for 0.4 < Y e < 0.5. 
A naive guess based on symmetry of NSE abundances 
would therefore not lead one to expect many free protons 
during freeze out for Y e > 0.5 either. A large number of 
free protons is however observed to occur in nuclear reac- 
tion networks calculations. We show that for a relatively 
cold, near freeze-out temperature NSE state, there is in- 
deed a qualitative difference in the mass fraction trends 
for Y e > 0.5 and Y e < 0.5. Restricting ourselves to 
the Ni isotopic chain and free nucleons, we explain how 
this difference arises as a result of a competition between 
the temperature dependent entropy term and the nuclear 
binding energy term contribution to the free energy. 

2. NUCLEAR STATISTICAL EQUILIBRIUM EQUATIONS 

NSE is established if all fusion reactions are in equi- 
librium with their inverses for a s et of thermody- 
nami c state variables and Y e (e.g. IClifford fc Tavlerl 
1965). Detailed balance relates the chemical po- 
tential of a nucleus Zi+ z Ni Xi to those of free nu- 
cleons: ^ = Zi/J,p + Ni/j, n . This yields an expres- 
sion for the number density of nucleus i: rii — 
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constraint equations can be perhaps most naturally writ- 
ten in a number density basis as yj 4 AiUi — ub and 
J2i Zm = n B Y ei where n B = n n + n p and Y e = n p /n n . 
We solve the NSE constraint equations numerically for 
the kinetic chemical potentials fi^ and /i^ m of the (as- 
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sumed Maxwellian) protons and neutrons respectively. 
For the Coulomb contributions /i^ and /if and the nu- 
clear partition func tions Qj we use the same form alism as 
ICalder et all (|2007l ) and (Seitenzahl et alj (|2008t) . 

3. RESULTS 

The results presented here are all for a baryonic mass 
density of p = 10 7 g cm~ 3 , which corresponds to a baryon 
number density of ub ~ 6.0 x 10 30 cm~ 3 . At high tem- 
perature (Tg = 9.0) the NSE mass fractions are domi- 
nated by free nucleons and 4 He. The 4 He mass fraction 
is symmetric across the line Y e = 0.5, and the mass frac- 
tions of free protons and neutrons are symmetric in a 
complementary sense - free protons are more abundant 
for Y e > 0.5 and free neutrons are more abundant for 
Y e < 0.5 (see fig.EQ). 

At somewhat lower temperature (Tg = 6.5), the sym- 
metry of the NSE mass fractions across the line of self- 
conjugacy is broken and only qualitatively discernible. 
For A Y" e > 0, free neutrons are less abundant for 
Y e — 0.5 — AY e than are free protons for Y e — 0.5 + AY e 
(see fig. [2]). Furthermore, the mass fraction of 4 He is 
not symmetric anymore and the abundance peaks of the 
Fe-peak nuclei are wider on the proton-rich side. 

At even lower temperature (Tg = 3.5), the qualita- 
tive features of the mass fractions of nuclei in NSE as a 
function of Y e at fixed density and temperature change 
dramatically in the transition from the neutron-rich to 
the proton-rich side. For Y e < 0.5, the mass frac- 
tion landscape is comprised of a sequence of overlap- 
ping abundance peaks (cf. e.g. IClifford fe Tavleij 119651 : 
lHartmann et~all 119851: iNadvozhin fe Yudinl l200l) . Fe- 
peak nuclei with a proton to nucleon ratio equal or close 
to the prescribed Y e of the ensemble and a large bind- 
ing energy per nucleon, q/A, are the most abundant nu- 
clei. For Y e > 0.5, the picture changes abruptly (see 
fig. [3]). The mass fraction distributions of the Fe-peak 
nuclei are no longer peaked, but rather either slowly ris- 
ing or falling. 56 Ni remains the most abundant nuclear 
species by mass all the way out past Y e = 0.6. The mass 
fraction of free protons continues to rise. 52 Fe shows a 
similar trend like 56 Ni, albeit at a smaller abundance 
level. There is only a slow rise in the mass fractions of 
proton-rich Fe-peak nuclei for increasing Y e , the abun- 
dance peaks for the Fe-peak nuclei with a proton to nu- 
cleon ratio equal to Y e that are so prominent on the 
neutron- rich side are absent. Various animations of NSE 
in proton-rich environments may be downloaded from 
http : / / cococubed . asu . edu/ code_pages/nse . shtml 

4. DISCUSSION 

In NSE, the Helmholtz free energy T = (U - Q) -TS 
is mi nimized with respect to th e nuclide mass fractions 
(e.g. INadvozhin fe Yudinl [2005). Before we discuss and 
compare the terms that make up T for different choices of 
compositions, it is instructive to review the dependence 
of q/A on the number of nucleons along an isotopic chain, 
the attributes of which ultimately are responsible for the 
lack of symmetry. 

4.1. The nuclear binding energy 

q/A on the proton-rich side for the Ni isotopic chain 
decreases rapidly towards proton drip, whereas it gen- 
tly increases to a maximum at A rnax = 62 before slowly 
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Fig. 1. — Mass fractions of nuclei in NSE as a func- 
tion of electron fraction for a constant density p = 10 7 g 
cm -3 (iij = 6.0 X 10 30 cm -3 ) and relatively high temperature 
T = 9.0 x 10 9 K. Shown are all nuclei with mass fractions larger 
than 10 -2 . The mass fractions exhibit large degree of (in the case 
of nucleons complementary) symmetry across Y e = 0.5. 
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Fig. 2. — Mass fractions of nuclei in NSE as a function 
of electron fraction for a constant density p = 10 7 g cm -3 
(n B = 6.0 X 10 30 cm" 3 ) and temperature T = 6.5 X 10 9 K. Shown 
are some abundant nuclei with mass fractions larger than 10~ . 
The symmetry across Y e = 0.5 is already broken. 
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Fig. 3. — Mass fractions of nuclei in NSE as a function 
of electron fraction for a constant density p = 10 7 g cm -3 
(n s = 6.0 X 10 30 cm -3 ) and temperature T = 3.5 X 10 9 K. Shown 
are some abundant nuclei with mass fractions larger than 10~ 5 . 
The mass fractions on either side of Y e = 0.5 exhibit qualitatively 
very different behavior. 
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Fig. 4. — q/A for Ni isotopes. Note the rapid decline on the 
proton-rich side of 56 Ni and the particularly large gap between 
te Ni and 55 Ni. 

falling off towards neutron drip (see fig. HJ. Other iso- 
topic chains, such as the one for Fe, look qualitatively 
very similar. The shape of the q/A-cmve can be quali- 
tatively understood by considering a simple liquid-drop 
mass formula taking the volume, surface, asymmetry, 
and Coulomb terms into account 

Q = a v - asA- 1 / 3 - a A (l - 2Z/A) 2 - a c Z 2 A^l z (3) 

For simplicity, we have neglected other terms such as 
pairing, Wigner, or residual interaction. The derivative 
with respect to nucleon number is given by 



S(A) = 



dQ 
dA 



= ^ A -^+^(l-2Z/A)+^ 



A 



A- 7 '\ 
(4) 



The second (asymmetry) term is positive on the proton- 
rich side, zero for 2Z = A, and negative on the neutron- 
rich side, which gives the g/A-curve its general concave 
up shape. It is, however, not entirely symmetric about 
the line 2Z = A. In fact, since 
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( 5 ) 

the absolute value of the slope due to the asymmetry 
term alone is larger on the proton-rich side. The first 
(surface) and third (Coulomb) term are always posi- 
tive, both monotonically increasing with A, which con- 
sequently shifts the most tightly bound isotope to the 
neutron-rich side and further increases the asymmetry 
in the slopes on either side of the maximum. This sim- 
ple argument shows the essence of how the surface and 
Coulomb terms in the nuclear mass formula leads to a 
more rapid fall off in q/A on the proton-rich side. 56 Ni is 
doubly magic and a Wigner N = Z nucleus, resulting in 
an increase in q/A compared to the simple mass formula 
above and the dominance of 56 Ni in cold, proton-rich 
NSE. 
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Fig. 5.— 27 56 Ni nuclei and 28 protons on the left, and 28 
55 Ni nuclei on the right. Both compositions have identical electron 
fraction Y e = 28/55. Even though pure 55 Ni has a marginally 
higher Q, the NSE state at p ~ 10 7 g cm -3 near freeze out is 
closer to the mix of 56 Ni nuclei and protons depicted schematically 
in the left panel. For almost equal Q, the composition with more 
particles and higher entropy is statistically preferred 




Fig. 6. — Differences in internal energy U, entropy contribution 
TS and binding energy Q between a mix of 56 Ni (subscript 56) 
and free protons (for Y e > 0.5) or neutrons (for Y e < 0.5) and a 
composition of pure Ni isotope (subscript *) for different Y e . 
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4.2. Free energy for different compositions 

To get a better understanding for the different low tem- 
perature behavior of abundance trends as a function of 
Y e for the proton and neutron-rich regimes, we consider 
a simplified system and restrict the composition to nu- 
clides from the Ni isotopic chain and free nucleons. For 
a concrete example, let us compare two compositions 
at Y e = 28/55 ~ 0.5091, one consisting of pure 55 Ni, 
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Fig. 7. — Helmholtz free energy T of compositions consisting 
of Ni isotopes and nucleons as a function of Y e . T = 3.5 X 10 9 K. 
and us = 6.0 X 10 30 cm~ 3 . For Y e > 0.5, a mix of protons and 
test isotope 56 Ni (thick blue line with dots) is favored. Other test 
isotopes (dotted lines) result in larger T . For Y e < 0.5, a pure 
isotope composition (thick red line with squares) has the lowest T . 
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and the other consisting of a mix of free protons and a 
"test" nucleus, here 56 Ni (see fig. The mean bind- 
ing energy per nucleon is given by Q = , where 
qi the binding energy of nucleus i and q n — q P = 0. 
Q of 28 55 Ni nuclei ( 28x4 %f 5 l McV = 8.497 MeV/nuc) is 
slightly larger than that of 27 56 Ni nuclei plus 28 pro- 
tons ( 2 27x56 + 28^i V = 8 - 486 MeV/nuc). NSE, however, 
favors the state with the lowest not the lowest Q. For 
Tg = 3.5 and ub = 6.0 x 10 30 cm -3 we know that pure 
55 Ni is disfavored and should have higher T . This is in- 
deed the AU between the two states is negligible 
and the entropy term TS is much larger for the state 
with the free protons (see fig. O, resulting in lower T . It 
is evident from fig. [5] that with 56 Ni as the test nucleus, 
the increase in the entropy term on the neutron-rich side 
for the mixed composition is more than compensated by 
the increase in binding energy by the pure state. On the 
other hand, for the proton-rich side the Q of both com- 
positions is nearly equal, and the higher entropy of the 
state with more particles makes the difference. 

Using th e data from the l atest atomic mass eval- 
uation of lAudi et al.l (120031) and an equatio n of 
state (iTimmes fc Arnettl ll999t iTimmes fc Swestvl 120001 
iFrvxell et al.ll2000l ). we compute and compare T for such 
compositions for the whole Ni isotopic chain, with each 
isotope as the test nucleus at a time. If Y e < Z/A of the 
test nucleus (in the example above and in fig. [6] this is 
56 Ni), we use free neutrons instead of free protons. 

For any given Y e , we can derive Q for a mix of a test 
nucleus A Ni with the appropriate number of free protons 
or neutrons. 



28 

if — > Y e , ( 



(6) 
(7) 



where g( A Ni) is the total binding energy of A Ni. T of all 
such compositions is shown as a function of Y e in fig. [71 
It clearly shows that for Y e < 0.5 a pure composition 
of the Ni isotope with Y e = 28 j A has the lowest J-, and 
that for Y e > 0.5 a composition consisting of free protons 
and 56 Ni minimizes T . Other choices for the test nucleus 
have larger T (dotted lines in fig. [7]) . 



5. CONCLUSIONS 

We have presented mass fraction trends for NSE in 
proton-rich environments. We have explained the, at first 
sight peculiar, high free proton mass fractions for Y e > 
0.5 and large mass fraction of 56 Ni even out to Y e = 0.6 
by considering a simplified model. Restricting ourselves 
to the Ni isotopic chain, we have explicitly shown that for 
Y e > 0.5 T is minimized by a state consisting of a mixture 
of free protons and the self-conjugate 56 Ni, whereas for 
Y e < 0.5 a state with a pure composition made up of 
the Ni isotope that has Z/A closest to the prescribed 
Y e is preferred. In reality, other isotopic chains are of 
course accessible. Other Fe-peak nuclides (especially the 
slightly neutron-rich even isotopes of Fe) compete with 
Ni isotopes for the most tightly bound nucleus with Z/A 
near Y e (see fig. EJ, giving the familiar NSE abundance 
pattern. 
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